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ABSTRACT 

Context. 

Aims. This study aims at setting observational constraints on the continuum and line core intensity contrast of network and faculae, 
specifically, their relationship with magnetic field and disc position. 

Methods. Full-disc magnetograms and intensity images by the Helioseismic and Magnetic Imager (HMI) on-board the Solar 
Dynamics Observatory (SDO) were employed. Bright magnetic features, representing network and faculae, were identified and the 
relationship between their intensity contrast at continuum and line core with magnetogram signal and heliocentric angle examined. 
Care was taken to minimize the inclusion of the magnetic canopy and straylight from sunspots and pores as network and faculae. 
Results. In line with earlier studies, network features, on a per unit magnetic flux basis, appeared brighter than facular features. 
Intensity contrasts in the continuum and line core differ considerably, most notably, they exhibit opposite centre-to-limb variations. 
We found this difference in behaviour to likely be due to the different mechanisms of the formation of the two spectral components. 
From a simple model based on bivariate polynomial fits to the measured contrasts we confirmed spectral line changes to be a sig- 
nificant driver of facular contribution to variation in solar irradiance. The discrepancy between the continuum contrast reported here 
and in the literature was shown to arise mainly from differences in spatial resolution and treatment of magnetic signals adjacent to 
sunspots and pores. 

Conclusions. HMI is a source of accurate contrasts and low-noise magnetograms covering the full solar disc. For irradiance studies 
it is important to consider not just the contribution from the continuum but also from the spectral lines. In order not to underestimate 
long-term variations in solar irradiance, irradiance models should take the greater contrast per unit magnetic flux associated with 
magnetic features with low magnetic flux into account. 

Key words. Sun: activity - Sun: faculae, plages - Sun: photosphere - Sun: surface magnetism 
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1. Introduction 

Photospheric magnetic activity is the dominant driver of vari- 
ation in solar irradiance on rotational and cyclical timescales 
( |Domingo et al. 2009} . Magnetic flux in the photosphere is partly 
confined to discrete concentrations of kilogauss strengths, gen- 
erally described in terms of flux tubes (Stenflo 1973; [Spruit &| 
|Roberts 1983 ). The brightness excess, or contrast relative to 
the Quiet Sun, of flux tubes is strongly modulated by their size 
and position on the solar disc (see Solanki 1993 , for a review). 
Within these magnetic concentrations, pressure balance dictates 
an evacuation of the interior and consequent depression of the 
optical depth unity surface ( Sprui t| 1 976| > . The horizontal extent 
influences the effect of radiative heating from the surroundings 
through the side walls on the temperature structure and contrast 
( |Spruit & Zwaan|1981[|Grossmann-Doerfh et al.|1994) . The po- 
sition on the solar disc changes the viewing geometry, and there- 
fore the degree to which the hot walls are visible and the appar- 
ent contrast (Steiner 2005 ). Models describing the counteracting 
effects on solar irradiance of dark sunspots, and bright network 
and faculae, characterizing the latter by the magnetic filling fac- 
tor (related to number density) and position have been successful 
in reproducing more than 90% of observed variation over multi- 
ple solar cycles fWenzler et al.|2006||Ball et al.|20"l2) . Other fac- 
tors, such as inclination, internal dynamics, phase of evolution 



and surrounding convective motions affect the brightness excess 
of a given flux tube, but become less important when consider- 
ing the overall behaviour of an ensemble as is the case with such 
models (Fligge et al. 2000). The same is assumed of flux tube 
size, which enters these models only very indirectly, though it is 
known to have a significant effect on contrast. 

Evidently, the robust reconstruction of solar irradiance vari- 
ation from models based on photospheric magnetic activity is 
contingent, amongst other factors, on a firm understanding of the 
radiant behaviour of magnetic elements, in particular the varia- 
tion with size and position on the solar disc. While the radiant 
behaviour of sunspots is relatively well known (Chapman et al. 
|1994| |Mathew et al.|2007]) a nd sufficiently described by current 
mode ls ( |Maltby et al.||1986[ |Collados et~aL|l994| |Unruh et af 
|1999| ), the converse is true of network and faculae, and consti- 
tutes one of the main uncertainties in current solar irradiance 
reconstructions. This is primarily due to the difficulty in observ- 
ing such small-scale features, the detailed structure of which are 
only starting to be resolved ( |Lites et al. 2004, Lagg et al. 2010) 
with instruments such as the Swedish 1-m Solar Telescope, SST 
(Scharmer et al. 2003) and the Im aging Magnetograph exper- 
iment, IMaX ([Martinez Pillet et al.||2011) on-board SUNRISE 
( |Solanki et al'||2010[ |Barthol et al'l|2011| ). As such, the rela- 
tionship between radiance and size cannot, as yet, be studied 
directly. It is however appropriate and more straightforward to 
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consider instead the relationship between apparent intensity con- 
trast and magnetogram signal. Apart from small-scale magnetic 



fields observed in the quiet Sun internetwork ( jKhomenko et al. 
|2003||Liteset al.|2008||Beck & Rezaei|2009| l, magnetic concen- 
trations carrying more than a minimum amount of flux exhibit 
similar field strengths regardless of size (Solanki & Schmidt 



1993t |Solanki et al.p999) . Flux tubes also tend towards verti- 
cal orientation due to magnetic buoyancy. As flux tubes exhibit 
a narrow range of magnetic field strengths and are largely verti- 
cal, on average the magnetogram signal at a given image pixel 
approximately scales with the proportion of the resolution ele- 
ment occupied by magnetic fields. Also, although the relation- 
ship between magnetogram signal and distribution of flux tube 
sizes is degenerate (a given magnetogram signal can, for exam- 
ple, correspond to either a single flux tube or a concentration of 
numerous smaller ones), flux tube size appears, on average, to 
be greater where the magnetogram signal is greater ( Orti z et ah] 
|20U2l 

Relatively few studies examining network and faculae con- 
trast variation with magnetogram signal and position on the so- 
lar disc have been reported in the literature. The majority of 
studies from the past two decades employed high-resolution 
(< 0.5 arcsec) scans made with ground-based telescopes. For 
example, the efforts of Topka et al. ( 1992 1997| l and Lawrence 



et~aT1 ( [T993]) with the Swedis h Vacuum Solar Telescope (SVST) 



and of Berger et al. (2007b with the SST These studies suf- 



fer from variable seeing effects introduced by the Earth's at- 
mosphere and poor representation of disc positions, a limita- 
tion imposed by the relatively narrow FOVs (field-of-view). 
Ortiz et al.](|2002|i and|Kobel et al.|(|2011|i repeated the work 



of |Topka et al.| ( |i992[ |1997| l and |Lawrence et all ( fT993] ) util- 
ising observations from spaceborne instruments, and in so do- 
ing avoided seeing effects. Ortiz et al. ( 2002| l employed full- 



disc continuum intensity images and longitudinal magnetograms 
from the Michelson Doppler Imager (MDI) on-board the Solar 
and Heliospheric Observatory (SoHO). While full-disc MDI ob- 
servations presented a more complete coverage of disc positions, 
allowing the authors to derive an empirical relationship relating 
contrast to heliocentric angle and magnetogram signal, the spa- 
tial resolution is significantly poorer than in the SVST studies (4 
arcsec versus > 0.3 arcsec). [Kobel et aL] ( |201 l| l examined the re- 
lationship between contrast and magnetogram signal near disc 
centre using spectropolarimetric scans from the Solar Optical 
Telescope (SOT) on-board Hinode (Kosug fet al.|2007) . In this 
instance, the spatial resolution (0.3 arcsec) is comparable. 

In this paper we discuss continuum and line core intensity 
contrast of network and facular elements from full-disc observa- 
tions by the Helioseismic and Magnetic Imager (HMI) on-board 
the Solar Dynamics Observatory (SDO) spacecraft ( Schou et al.| 



2012 ), and their relationship with heliocentric angle and magne- 



togram signal. The aim here is to derive stringent observational 
constraints on the relationship between intensity contrast, and 
position on the solar disc and magnetic field. This will be of util- 
ity to solar irradiance reconstructions, especially as HMI data 
will increasingly be used for this purpose. 

This study partly echoes the similar studies discussed above, 
in particular that by |Ortiz et aL ([2002 ) utilising MDI observa- 
tions. It presents a significant extension of the effort by Ortiz 



et al. ( 2002 1 in that we examined the entire solar disc not just 



in the continuum but also in the core of the HMI spectral line 
(Fe I 6173 A). This is of particular relevance to solar irradiance 
variation given the observation that spectral line changes appear 
to have a significant influence on such variations ([M itchell & 
|Livingston|199T]|Unruh et al.|1999||Preminger et al.|2002) . Both 



here and in the study by |Ortiz et aL] ( |2002| l, network and facu- 
lar elements were distinguished from quiet Sun by the magne- 
togram signal, and sunspots and pores by the continuum inten- 
sity. HMI magnetograms are significantly less noisy than MDI 
magnetograms, allowing us to achieve a similar magnetogram 
signal threshold while averaging over a much shorter period 
(315 seconds versus 20 minutes). Network and facular features 
evolve at granular timescales (~ 10 minutes, Berger & Title 



1996; Wiegelmann et al. 20121. It is pertinent to keep the av- 
eraging period below this in order to avoid smearing and loss 
of signal. HMI also has a finer spatial resolution (1 arcsec com- 
pared to 4 arcsec), allowing weaker unresolved features to be 
detected at the same noise level than with MDI. The finer res- 
olution however, also renders intensity fluctuations from small- 
scale phenomena such as granulation and filamentation more se- 
vere, which complicates the clear segmentation of sunspots and 
pores. 

In Sect. |2.1| we briefly present the HMI instrument, the ob- 
servables considered and the data set. The data reduction process 
by which we identified and derived the intensity contrast of net- 
work and facular features is detailed in Sect. |2.2| Following that 
we describe the results of our analysis of these measured con- 
trasts (Sect.[3]l. In Sect. |4] we discuss our findings in the context 
of earlier studies and of their relevance to facular contribution to 
solar irradiance variation, before presenting our conclusions in 
Sect. [5] 

2. Method 

2.1. SDO/HMIdata 



SDO/HMI (Schou et al. 2012 1 is designed for the continuous, 
full-disc observation of velocity, magnetic field and intensity on 
the solar surface. The instrument comprises two 4096 x 4096 
pixel CCD cameras observing the Sun at a spatial resolution of 
1 arcsec (corresponding to two pixels). By means of a tunable 
Lyot filter and two tunable Michelson interferometers, the in- 
strument records 3.75-second cadence filtergrams at various po- 
larizations and wavelengths across the Fe I absorption line at 
6173 A. 45-second cadence Dopplergrams, longitudinal magne- 
tograms and intensity (continuum, line depth and width) images 
are generated from the filtergram sequence. For this work we 
considered the longitudinal magnetic field, continuum intensity 
and line depth observables. HMI is full-Stokes capable, however, 
at time of study, only 720-second cadence Stokes IQUV param- 
eters and Milne-Eddington inversions were available. As argued 
in Sect. [T] for this study it is important to keep the integration 
period of measurements below ~ 10 minutes. We opted to utilise 
the 45-second longitudinal magnetograms also to keep in line 
with earlier studies, which examined intensity contrast variation 
with line-of-sight magnetic field (e.g.|Topka et al.|[T992, 1997 



|Lawrence et al.|1993||Ortiz et al.|2002||Kobel et al.|201 1). M ore 
details on the instrument can be found in |Schou et al.| ( |2012"[ ). 

The data set comprises simultaneous longitudinal magne- 
tograms, continuum intensity and line depth images from 15 
high activity days in the period May 2010 to July 2011. From 
each day, for each observable we took the average of the seven 
45-second cadence images from a 315-second period, each ro- 
tated to the observation time of the middle image to co-register. 
Aside from signal-to-noise considerations, this averaging is to 
suppress variance from p-mode oscillations. The dates and times 
of the employed observations are listed in Table [T] 

Longitudinal magnetograms describe the line-of-sight com- 
ponent of the average magnetic flux density over each resolution 
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Table 1. Observation date and time of the data set. 



Date 


Time 


(year.month.date) 


(hour.minute. second) 


2010.05.04 


00:03:00 


2010.06.11 


00:00:00 


2010.07.24 


00:05:15 


2010.08.11 


00:00:00 


2010.09.02 


00:00:00 


2010.10.25 


00:00:00 


2010.11.13 


00:00:00 


2010. 12.04 


ria, AA. AA 

00:00:00 


2011.01.01 


00:00:00 


2011.02.14 


00:00:00 


2011.03.08 


00:00:00 


2011.04.15 


00:00:00 


2011.05.30 


00:00:00 


2011.06.02 


00:00:00 


2011.07.18 


00:12:00 



Notes. The time listed is the nominal time, in International Atomic Time 
(TAI), of the middle cadence in the sequence of seven considered from 
each day. 



element. To first order, for a given flux tube of intrinsic magnetic 
field strength B, the unsigned longitudinal magnetogram signal, 
(B[), is aB\ cos y\, scaled by a, the magnetic filling factor and y, 
the inclination of the magnetic field from the line-of-sight. As 
mentioned in the introduction (Sect. fTJ, flux tubes tend towards 
vertical orientation. In this study we examined the overall prop- 
erties of an ensemble of magnetic elements. Under this condition 
it is reasonable to assume that on average y m 0, the heliocentric 
angle, allowing us to employ the magnetogram signal as a proxy 
for aB via the quantity (B/) /fi, where /i = cosfl. This quantity 
also represents a first-order approximation of the unsigned aver- 
age magnetic flux density over each resolution element. At time 
of writing, HMI Dopplergrams and longitudinal magnetograms 
were generated from just the first Fourier component of the filter- 
gram sequence resulting in a V2 factor increase in photon noise 
from the optimal level ( |Couvidat|2011 Couvidat et al.|2012) . 

Each line depth image was subtracted from the correspond- 
ing continuum intensity image to yield the line core intensity im- 
age. Hereafter we will denote the continuum and line core inten- 
sity Ic and /l respectively. At time of writing, HMI data products 
are generated from the filtergram sequence assuming a Gaussian 
form to the Fe I 6173 A line and delta filter transmission pro- 
files. The effect of these approximations on Doppler velocity 
and longitudinal magnetic field measurements are accounted for, 
but not completely for the intensity observables ( |Couvidat]2 011 ; 
Couvidat et al.|2012 1. The impact on this study is assumed min- 
imal as we are only interested in contrast relative to the local 
mean quiet Sun level. 



2.2. Data reduction 

2.2.1. Magnetogram noise level 

The noise level of 315-second HMI magnetograms as a function 
of position on the solar disc was determined. For this purpose 
we used 10 spot-free 315-second magnetograms recorded over 
a seven month period in 2010. First, we estimated the centre-to- 
limb variation, CLV of the noise level. The pixels within each 
magnetogram were ordered by distance from the disc centre and 
sampled in successive blocks of 5000 pixels. The blocks repre- 
sent concentric rings (near the limb, arcs, as the circumference 



of the solar disc is greater than 5000 pixels) of pixels of similar 
distance from the disc centre. Within each block we computed 
the average distance from the disc centre and the standard devi- 
ation of the magnetogram signal (iteratively, with points outside 
three standard deviations from the mean excluded from the suc- 
ceeding iteration till convergence). A fifth-order polynomial in 
fi was fitted to the noise CLV; the mean of the standard devia- 
tion versus distance profiles so derived from the magnetograms. 
The magnetograms were then normalized by the noise CLV fit. 
At each disc position, the standard deviation over a 401 x 401 
pixel window centred on the pixel of interest was computed (it- 
eratively as above) for each normalized magnetogram and the 
median taken (following |Ortiz et al. 2002). (Near the limb, the 
standard deviation was computed from just the image pixels that 
lie within the solar disc.) A sixth-order polynomial was fitted to 
the resultant surface (termed here the noise residue). This fit rep- 
resents the noise level after the removal of the CLV as a function 
of position on the solar disc. The noise level, cr^,), shown in Fig. 
[T] is then the product of the noise CLV fit and the noise residue 
fit. The noise level is lowest near disc centre and increases radi- 
ally up to the limb (mean of 4.9 G for fi > 0.95 and 8.6 G for 
fi < 0.05). The root-mean-square, RMS difference between the 
noise level and the noise CLV fit is 0.4 G. The correlation with 
distance from disc centre and relatively small deviation from cir- 
cular symmetry suggests photon noise is the dominant compo- 
nent. The noise level of 45-second magnetograms, derived by a 
like analysis, has a similar, albeit accentuated form. The ratio be- 
tween the noise level, averaged over the solar disc, of 45-second 
and 315-second magnetograms is 2.7 (approximately V7). The 
noise level of 45-second magnetograms determined by Liu et al. 
(2012), via a vastly different method, also exhibits a similar CLV. 
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Fig. 1. Magnetogram noise level in units of Gauss, as a function 
of disc position, sampled at 16-pixel intervals in either direction 
on the CCD. 



2.2.2. Quiet Sun intensity 

For this part of the data reduction process, where we examined 
how quiet Sun intensity and the noise level of the intensity im- 
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ages vary with position on the solar disc, we counted all pixels 
with (B/) /fi < 10 G as corresponding to quiet Sun (QS). 

The continuum and line core intensity images were normal- 
ized by the fifth-order polynomial in // fit to the quiet Sun pixels 
to correct for limb-darkening (following Neckel & Labs 1994). 
In the case of the line core intensity images, this was also to 
correct for the centre-to-limb broadening of the Fe I 6173 A 
line (Norton et al.|[2006). There are distortions in the intensity 



images such that after this normalization, the mean quiet Sun 
intensity is not constant at unity but varying with position on 
the solar disc. The mean quiet Sun intensity at the continuum 
and line core as a function of position on the solar disc, denoted 
(7c,qs) and (^l,qs) respectively, were determined for each of the 
selected days. A 401 x 401 pixel window was centred on each 
disc position and the mean continuum and line core intensity of 
all quiet Sun pixels inside the window computed. For each day 
represented in the data set, (/c,qs) and (/l,qs) were given by the 
fifth-order polynomial fits to the mean quiet Sun continuum and 
line core intensity surfaces so derived from the images from the 
day. This analysis had to be repeated for each selected day as 
we found the spatial distribution of (/c,qs) and (/l,qs) to vary 
significantly over the period of observation. The RMS value of 
(Ic.qs) - 1 and (/l,qs) - 1. the scale of the image distortions, is 
on average 0.004 and 0.01 respectively. 

The CLV of the standard deviation of quiet Sun intensity at 
continuum, cr IcQS and line core, ct /lqs were determined, from 
Icl (Ic,qs) and 1^/ (Il,qs), by an analysis similar to the proce- 
dure used to derive the magnetogram noise CLV. In Fig. [2j we 
express cr /cQS and o"/ LQS , which carry information on the noise 
level of the intensity images and granulation contrast, as a func- 
tion of fi. Going from disc centre, <x /cQS decreases gradually 
down to ju ~ 0.15 before increasing rapidly towards the limb. 
The monotonic decline from disc centre to fi ~ 0.15 resembles a 
similar trend in the CLV of granulation contrast reported by var- 
ious authors (Sanchez Cuberes et al. 2000 2003 and references 
therein). ct /lqs exhibits a similar, though less accentuated, trend. 
For both cr /cQS and ct /lqs , the /v-dependence from disc centre to 
/i ~ 0.3 is approximately linear, as highlighted by the linear fits 
(dotted lines). The elevation near limb is a direct consequence of 
limb-darkening; the diminishing signal-to-noise ratio translates 
into an escalating noise level in the normalized intensity. Given 
the gross scatter towards the limb, we excluded pixels outside 
jj = 0.1, representing about 1% of the solar disc by area, from 
further consideration. 



2.2.3. Identification of network and faculae 

Network and facular features, the subject of this work, were 
identified by first distinguishing them from quiet Sun by the 
magnetogram signal and from sunspots and pores by the con- 
tinuum intensity. Pixels with {B/} In > 30-^/fj. (~ 14 G near 
disc centre, where o-(b,) is lowest) were considered to harbour 
substantive magnetic fields. Isolated pixels above this threshold 
were assumed false positives and excluded. Hereafter we will de- 
note 3<T(b,)I^ as ((Bi) //i) TH - Pixels with Iq < 0.89 were counted 
as sunspots and pores. This value of the threshold is given by 
the mean of the minimum value of (/c,qs) _ 3cr /cQS from each 
selected day. It was so defined to distinguish sunspots and pores 
with minimal wrongful inclusion of intergranular lanes and mag- 
netic features darker than the quiet Sun in the continuum. The 
continuum intensity versus yu scatter plot of network and facular 
pixels for one of the selected days (June 2, 2011) is shown in 
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Fig. 2. Standard deviation of quiet Sun ((B/> /fi < 10G) intensity 
at continuum, cr /cQS (red) and line core, ct /lqs (blue) as a function 
of \i. The dotted lines represents the linear fit to cr /cQS and 07 
over the range 0.3 < fi < 1.0 (largely hidden due to the close 
agreement), extrapolated to yu = 0. The dashed line denotes the 
threshold (yu = 0.1) below which pixels were excluded from the 
rest of the study in view of the scatter in measured intensity. 



Fig. [3] (top panel). The pixels identified as network and faculae 
lie clearly above the continuum intensity threshold. 

In HMI magnetograms, magnetic signals produced by 
sunspots and pores extend beyond their boundary (in our analy- 
sis, the Ic = 0.89 locus). This is illustrated for a sunspot near disc 
centre and another that is close to the limb on one of the selected 
days (July 18, 2011) in Fig. [4] The I c = 0.89 locus is plotted 
over the continuum intensity image and magnetogram of both 
sunspots (red contours) to highlight the extension of the magne- 
togram signal. This arises predominantly from the lateral expan- 
sion of their magnetic field with height (i.e., magnetic canopy) 
and partly also from the effect of straylight from sunspots and 
pores on nearby pixels. Towards the limb, these magnetic fea- 
tures become more extensive and bipolar due to the acute orien- 
tation of the largely horizontal magnetic canopy with the line-of- 
sight ( Giovanelli|198 0), as illustrated by the near limb example 
(bottom panels). Figure [4] also highlights the presence of bright 
sunspot structures (such as bright penumbral filaments) that lie 
above the continuum intensity threshold. All of these could eas- 
ily be misidentified as network and faculae by the simple ap- 
plication of the magnetogram signal and continuum intensity 
thresholds described earlier. 

To account for the effects of straylight around sunspots and 
pores, magnetic canopies and bright sunspot structures, we ex- 
panded sunspots and pores to include adjoining magnetic signal. 
Pixels adjacent to sunspots and pores that lie above ((Bi) /yu) TH 
are reassigned to these features. This was iterated till no more 
pixels could be added. Here we will refer to this process as mag- 
netic extension removal. Simply adding only adjoining pixels 
within a threshold distance from sunspots and pores instead is 
not useful as the physical extent of magnetic canopies exhibits 
a broad dynamic range, dependent on the position and physi- 
cal properties of the associated feature. Kobel et al. (201 1|>, in a 



similar study with Hinode/SOT scans, expanded pores to include 
adjoining pixels above a threshold magnetogram signal level of 
200 G to account for their influence on surrounding pixels from 
telescope diffraction. While this technique appears to work for 
Hinode/SOT scans, here we observed that adding only adjoin- 
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Fig. 3. Continuum intensity, Iq versus fi (top) and (Bi) /fi versus 
jj (bottom) scatter plots of pixels counted as network and faculae 
(green) from June 2, 2011. The pixels counted as sunspots and 
pores by the continuum intensity threshold and the magnetic ex- 
tension removal procedure (see text) are represented by the cyan 
and blue dots respectively. The red dots represent network and 
faculae pixels that lie above the cutoff {Bj) /fi level. The dashed 
lines denote the continuum intensity threshold (top) and the cut- 
off (Bi) In level (bottom). 



ing pixels above a threshold magnetogram signal level results in 
the appearance of a knee in contrast versus {Bf) /fi plots at this 
threshold level regardless of the value chosen. For these reasons 
we expanded sunspots and pores to include all adjoining mag- 
netic signal, though this conservative approach inevitably as- 
signs too many pixels, including legitimate faculae, to sunspots 
and pores. 

All pixels identified as magnetic, and not as sunspots and 
pores or their extensions, were taken to correspond to network 
and faculae. Figure [5] shows the continuum intensity versus fi 
(top panel) and (B/) /fi versus fi (bottom panel) scatter plots of 
pixels identified as network and faculae (green), and counted as 
sunspots and pores by the magnetic extension removal procedure 
(blue) from June 2, 2011. The pixels captured by the magnetic 
extension removal procedure are not well distinguished from 
network and faculae by the continuum intensity; largely hidden 
by network and faculae in the continuum intensity versus fi scat- 
ter plot. It is clear from the {Bi) /fi versus /j. scatter plots however 
that the two classes are significantly different magnetic popula- 
tions. As noted earlier, this procedure is likely too severe and 
results in the exclusion of some true faculae. However, for the 
purpose of this study it is not necessary to identify all faculae 
present and far more important to avoid false positives. 

Finally, network and facular pixels with (Bi) /fi above a con- 
servative cutoff level of 800 G were excluded from the sub- 
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Fig. 4. Continuum intensity image (top left) and magnetogram 
(top right) from July 18, 2011 and the corresponding 50 x 
50 arcsec insets of the boxed sunspot features near disc centre 
(middle panels) and limb (bottom panels). The red contours rep- 
resent the continuum intensity sunspot boundary (the Iq = 0.89 
locus). The magnetograms have been saturated at +80 G to high- 
light the extension of the magnetic signal from these sunspots 
beyond the intensity boundary (see grey scale below the figure). 



sequent analysis (following Ball et al. 201 1| [20T2 1. They are 
mostly bright features concentrated near the limb (as illustrated 
for June 2, 2011 by the red dots in Fig. associated with 
sunspots and pores. (The relatively high {Bf) jyL values likely re- 
flect nearly horizontal fields, for which | cos y\ » fi towards the 
limb.) This is to account for non-facular magnetic signals that 
might have been missed by the continuum intensity threshold 
and the magnetic extension removal procedure. 

The classification image, indicating the positions of the pix- 
els classed as sunspot and pores, and network and faculae for 
another of the selected days (May 30, 2011) is shown in Fig. [5] 
In spite of the severe measures taken to minimise the influence 
of sunspots and pores, a fair fraction of active region faculae re- 
mains. In total, 7.6 x 10 5 pixels were identified as corresponding 
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to network and facular features from the data set (i.e. 4.5% of all 
solar disc pixels in the 15 continuum intensity image and mag- 
netogram pairs examined). 

2.2.4. Definition of intensity contrast 

The average continuum and line core intensity contrast over a 
given resolution element, C/ c and C/ L were defined as: 



C/ c (x,y) = 



and 



C/ L ( x >y) = 



Ic(x,y) - (I c ,Qs)(x,y) 
(Ic,Qs)(x,y) 



I L (x,y) - </ L ,Qs)fey) 
(Ih,Qs)(x,y) 



(1) 



(2) 



respectively, where (x,y) denote position on the CCD array. 
These two quantities, computed for each of the pixels identified 
as corresponding to network and facular features, represent the 
normalized difference between the continuum and line core in- 
tensities at a given pixel and the local mean quiet Sun levels as 
given by the mean quiet Sun continuum and line core intensity 



surfaces, (/c,qs) an d (^l.qsX defined earlier in Sect. 2.2.2 



In summary, here we extracted an ensemble of 7.6 x 10 6 
continuum and line core intensity contrast measurements cor- 
responding to network and facular features covering as wide a 
range of heliocentric angles (0.1 < n < 1.0) and magnetogram 
signal (((-S/) //0th < (®z) /M — 800 G) as reasonably possible 
from the data set for the succeeding analysis. 



3. Results 

3.1. Variation with position and magnetogram signal 

The positions of pixels identified as network and faculae, clas- 
sified by (B/) In, in a quiet region and an active region on one 
of the selected days (April 15, 2011) is shown in Fig [6] As ex- 
pected, at HMI's spatial resolution, magnetic signals with higher 
{Bi) //j are largely concentrated in active regions. Though mag- 
netic signals with lower (B/) ///are present in both quiet Sun and 
active regions, the fact that the solar disc is predominantly quiet 
Sun means these signals correspond largely to quiet Sun network 
and internetwork. 

To elucidate the CLV of intensity contrast, we grouped the 
measured contrasts into eight intervals of {Bi) In spanning the 
range ((B/) /aOth < (fii) /A* - 800 G and within each inter- 
val into ju bins 0.05 wide. As the distribution of magnetogram 
signal is skewed towards the lower bound ( Wenzler et al.| 2004 



Parn ell et al.|2009] l, the {Bi) In intervals were defined such that 
the widths slide from about 36 G {({Bi) ///) TH < {B,} In < 50 G) 
to 160G (640G < (B t ) In < 800G) to ensure reasonable statistics 
in every interval. In grouping the measured contrasts into these 
broad {B/} In intervals we are effectively grouping the network 
and facular features by aB (Sect. 2.1 1, neglecting differences in 
quiet Sun network and active region faculae contrast ( [Lawrence] 
|et al.| 1993 1 Ko bel et al.|201 \\ . This is reasonable since the lower 
intervals are mainly populated by quiet Sun network and the 
higher intervals by active region faculae. The bin-averaged con- 
tinuum and line core intensity contrast as a function of n and 
the cubic polynomial fit for each of the {Bi) In intervals are ex- 
pressed in Figs.|7]and|8]respectively. For brevity we will refer to 
these bin-averaged contrasts as the contrast CLV profiles. 

The truncated /i coverage in the lower {Bfi //i intervals is due 
to foreshortening. As flux tubes are mainly vertical, going from 



disc centre to the limb the corresponding longitudinal magne- 
togram signal diminishes, eventually dropping below the thresh- 
old level which itself rises towards the limb (Fig. [TJ. The fluc- 
tuations near the limb, more pronounced in the higher {Bi) In 
intervals, is due to the inhomogeneous distribution of active re- 
gions on the solar disc, where stronger magnetic signals are con- 
centrated on the selected days. Diminishing statistics also play a 
role; there are comparatively fewer pixels in the higher {Bf) In 
intervals. 

To investigate the {Bi) //i-dependence of intensity contrast, 
the measured contrasts were grouped into eight n intervals span- 
ning the range 0.1 < n < 1.0 and within each interval into {B\) In 
bins 40 G wide. The n intervals were defined such that they rep- 
resent an approximately equal proportion of the solar disc by 
area. The bin-averaged continuum and line core intensity con- 
trasts as a function of {Bi) In and the cubic polynomial fit for 
each of the n intervals are shown in Figs. [9] and 10 respectively. 
For brevity we will refer to these bin-averaged contrasts as the 
contrast versus {B{) In profiles. These profiles represent the vari- 
ation of continuum and line core intensity contrast ranging from 
internetwork and weak network to active region faculae at dif- 
ferent distances from disc centre. 

The cubic polynomial fit to each contrast versus {Bi) In pro- 
file included the zeroth-degree term. This produced better fits to 
data at low {Bi) In than constraining the fits to pass through the 
origin by excluding it. Approaching {Bj) In — G, the contrast 
versus {Bi) In profiles decline to weak, broadly negative levels. 
Although the cubic polynomial fits express good agreement with 
measurements, they are too simple to accommodate this decline 
well without including the unphysical zeroth-degree term. For 
similar considerations we included the zeroth-degree term in the 
cubic polynomial fits to the contrast CLV profiles, so rendering 
them non-zero at n — 0. These apparent offsets in the contrast 
CLV and contrast versus {Bi) In profiles probably reflects the 
fact that magnetic elements are generally located in dark inter- 
granular lanes. 

In Fig.[TT]we show a recomputation of the continuum inten- 
sity contrast versus {Bi) In profile over 0.94 < < 1.00 from 
Fig. 9] where we included pixels below the magnetogram sig- 
nal threshold and not identified as sunspots and pores (i.e., quiet 
Sun), and grouped the measurements by {Bj) In in bins 10 G 
(instead of 40 G) wide. Approaching {Bi) In = G, contrast de- 
clines gradually to below the reference level before turning back 
up sharply towards the origin. Similar trends were reported by 
I Narayan & Scharmer| \20\0) , |Kobel et al] ( |20 1 1 [ > and |Schnerr 
|& Spruit ( |2011| i, who termed it the fishhook feature, based on 
SST and Hinode/SOT disc centre scans. Schnerr & Spruit ( 201 \) 
demonstrated the resolution of granules and dark intergranular 
lanes, where magnetic flux concentrates, at the relatively fine 
spatial resolution of both instruments (0.15 and 0.3 arcsec re- 
spectively) to be the cause of the fishhook feature at low mag- 
netogram signal levels. Though HMI has a coarser resolution (1 
arcsec) than either Hinode/SOT or SST, the fishhook feature near 
{Bi) In — G in Fig. 1 1 indicates granulation is still sufficiently 



resolved to have a measurable impact on apparent contrast. 

Comparing the measured intensity contrast at continuum and 
line core derived here, the most distinct difference is the opposite 
CLV. Confining this discussion to broad trends in the measure- 
ments, continuum intensity contrast is weakest near disc centre 
and increases up to a maximum before declining quite signifi- 
cantly towards the limb (Fig.[7]l. Conversely, line core intensity 
contrast is strongest near disc centre and declines monotonically 
from disc centre to limb (Fig. |HJ. Line core intensity is mod- 
ulated by line strength and shape, and continuum intensity. The 
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Fig. 5. HMI continuum intensity image (left), magnetogram (middle) and classification image (right) from May 30, 2011. The 
classification image indicates the positions of the pixels classed as network and faculae (white), and sunspots and pores (black). 
The latter includes magnetic signal adjoined to sunspots and pores, counted to them to avoid counting their magnetic canopy, 
possible bright structures within penumbrae and straylight as network and faculae erroneously. The continuum intensity image and 
magnetogram are scaled between 0.6 and 1.2, and -80 G and 80 G respectively as in Fig. [4] 
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Fig. 11. Continuum intensity contrast, C/ c of quiet Sun, network 
and faculae over 0.94 < /j. < 1.00 as a function of {Bi) /fi. The 
filled circles represent the mean of C/ c grouped by {Bi) ///in bins 
10 G wide. The dashed line denotes the mean quiet Sun level. 



centre-to-limb decline exhibited by the measured line core inten- 
sity contrast arises from variation in line strength and shape (if 
this were absent the continuum and line core intensity contrast 
would vary proportionally) and would be even more acute than 
reflected in the CLV profiles if not partially offset by the accom- 
panying increase in continuum intensity. This will be demon- 
strated, along with a closer discussion of the diverging trends 



exhibited by both sets of measurements in Sect 4.3 



3.2. Surface fits 

The cubic polynomial fit to the contrast CLV and contrast versus 
(B/) /fi profiles reproduced the observations well (Figs. [7]to[T0]i. 
Given this, we fit the entire set of measured network and fac- 
ulae continuum and line core intensity contrast as functions of 
fi and {Bi) jfi (following |Ortiz et~a l. 2002). The measured con- 
trasts were grouped by fi and the natural logarithm of {Bi) jfi 
into a grid of 36 x 41 bins. The grid spans 0.1 < /J < 1.0 and 
2.6 < ln«fi/> /ft) < 6.7 or about 14 G to 810 G. Each bin repre- 



sents an interval of 0.025 in fi and 0.1 in ln((B/) /fi). The loga- 
rithmic binning in (Bi) //i is to compensate for the bottom-heavy 
distribution of magnetogram signal, to ensure the distribution of 
points is not too concentrated in the lower magnetogram signal 
bins. At each grid element we considered the mean fi, (B/) ///, 
continuum and line core intensity contrast of the points within 
the bin. In accord with the cubic polynomial fits to the individual 
contrast CLV and contrast versus (Bi) /fi profiles, bivariate poly- 
nomials cubic in fi and (B/) /fi were fitted to the bin-averaged 
continuum and line core intensity contrast. The zeroth fi and 
{B{) j fi orders were included based on similar considerations as 
with the polynomial fits to the contrast CLV and contrast versus 
(Bi) jfi profiles. In linear algebra notation, the surface fit to the 
bin-averaged continuum and line core intensity contrast are: 







10 




10 




10 




10 



-5.11 
-0.04 
0.19 
-0.08 



7.74 
3.84 
-6.27 
3.58 



0.34 
-7.42 
12.03 
-8.04 



-4.72 
3.90 

-6.78 
5.04 



(3) 



and 



C 



io- 
io- 



10 -8 (<|>y 



-1.08 
-0.08 
0.07 
-0.06 



2.09 
6.84 
-1.50 
1.03 



-1.78 
-11.22 

2.52 
-1.90 



0.66 
6.34 
-1.48 
1.18 



(4) 



respectively. Since contrast is wavelength dependent and bright 
magnetic features are largely unresolved at HMI's spatial reso- 
lution, these relationships are valid only at the instrument's op- 
erating wavelength (6173 A) and spatial resolution (1 arcsec). 

The surface fits are illustrated as surface and grey scale plots 
in Figs. [T2"| and [T3| Cross-sections to the surface fits are plotted 
in Figs. |7|to[T0| (blue curves) along the contrast CLV and contrast 
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versus (Bj) //x profiles, and the corresponding cubic polynomial 
fits (red curves). The surface fits are in excellent agreement with 
the cubic polynomial fits almost everywhere. The agreement is 
so close that the surface fit cross-sections are completely hidden 
by the cubic polynomial fits in most places. The advantage with 
these surface fits is that they allow us to describe how the mea- 
sured contrasts vary with fi and (B{) Ifi almost equally well and 
with far fewer free parameters than by the cubic polynomial fit to 
each individual contrast CLV and contrast versus {Bi) Ifi profile 
(32 versus 128 free parameters; 4 from each of 32 profiles). 
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Fig. 12. Grey scale plot of the polynomial fit to bin-averaged 
continuum intensity contrast over the region where data exist, 
sampled at the mid-point of each bin in the 36 x 41 bins grid 
employed to compute the averages. 
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Fig. 13. Same as Fig. 12 but for line core intensity contrast. 



3.3. Intrinsic contrast 

At each (Bf) Ifi interval, we estimated the maximum continuum 
intensity contrast, C/ Cjlnax and the heliocentric angle at which 
it is reached, fi lmx from the cubic polynomial fit to the con- 
trast CLV profile (Fig. [7). Following |Ortiz et al.] ( |2002} we 
term C/ Cjlnax / ((Bf) Ifi) the specific contrast; estimated here from 
the quotient of C/ C . max and ((Bi) Ifi), the mean (B/) Ifi of all 



the points within a given (Bi) Ifi interval. The values of yu max , 
C/c.max and C/ C>max / ((Bf) Ifi) are shown in Fig. fl4]as a function 
of (Bi) Ifi where the abscissa is given by ((Bf) Jfi). Also plotted 
are the values obtained from the surface fit (Equation|3]l. The un- 
certainty in C/ C , m ax is given by the RMS difference between the 
contrast CLV profiles and their cubic polynomial fits. The uncer- 
tainty in C/ Cimax / ((Bi) Ifi) was estimated from the uncertainty in 
C/ C , max and the standard error of ((B/) Ifi) employing standard 
propagation of errors. 

While the position of the continuum intensity contrast CLV 
maximum varies with (B/) Ifi, line core intensity contrast at 
a given (B/) Ifi is invariably strongest at disc centre (Fig. [8]l 
as pointed out in Sect. 3.1 Maximum line core intensity con- 
trast (i.e., the value at = 1), C/ Li/J= i and specific contrast, 
C/ L ,/i=i / ((Bi) Ifi) derived similarly as above from the cubic poly- 



nomial fit to the contrast CLV profiles are expressed in Fig. 15 
as a function of (Bf) Ifi, along with the values obtained from the 
surface fit (Eq. [4]). In both instances, there is some disparity be- 
tween the values obtained from the cubic polynomial fits to the 
contrast CLV profiles and the surface fit below (Bf) jfi ~ 100 G. 
This is likely due to the truncated coverage of disc positions at 
low (B/) Ifi discussed in Sect. 3.1 
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Fig. 14. Heliocentric angle at which continuum intensity contrast 
reaches its maximum, // max (top), as well as the contrast, C/ Cimax 
(middle) and specific contrast, C/ Cjmax / ((Bi) In) (bottom) there, 
as a function of (B{) jfi. The filled circles represent the values 
derived from the cubic polynomial fit to the contrast CLV pro- 
files (Fig. [7| and the error bars the uncertainty in C/ C>max and 
C/ Cjlnax / ((Bf) Ifi). The curves follow the solution from the sur- 
face fit to measured continuum intensity contrast (Eq.[3]l. 

Network and facular features exhibit similar, kilogauss 
strength magnetic fields ( |Stenflo|1973||Solanki & Stenflo|1984l 
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Fig. 15. Line core intensity contrast, C/ L „ = i (top) and specific 
contrast, C/ LiM= i/ ((Z?/) ///) (bottom) at disc centre as a function 
of (Bi) fix. The filled circles represent the values derived from the 
cubic polynomial fit to the contrast CLV profiles (Fig. [8} and the 
error bars the uncertainty in C/ L ^=i and C/ Li/J= i/ ((Z?/) ///). The 
curves follow the solution from the surface fit to measured line 
core intensity contrast (Eq.Q. 



Rabi n|1992| Riiedi et al.|1992 1. This means, in general, the in- 
trinsic magnetic field strength of magnetic elements, B does not 
vary significantly with (Z?/) //i and is also well above the range of 
{Bi) /// examined here, thus allowing us to assume that the mag- 
netic filling factor, a never saturates. As mentioned in Sect. 2.1 
the quantity {Bi) /fi is representative of aB. Taking these into 
account, within the limits of this study, {Bi) jfi is representative 
of the magnetic filling factor. Therefore specific contrast, as it 
is defined and derived here, is a measure of the average intrin- 
sic contrast or heating efficiency of magnetic features within a 
given {B{) jfi interval. Like the intrinsic field strength, the intrin- 
sic contrast cannot be measured directly because most network 
and facular features are not resolved at HMFs spatial resolution 
and due to the effect of scattered light. 

Above {Bi) //J ~ 100 G, the heliocentric angle at which con- 
tinuum intensity contrast reaches its maximum decreases (i.e., 
the maximum is reached further away from disc centre) with in- 
creasing (Bi) In (Fig. 14 1. A similar trend was reported by Ortiz 
et al. ( 2002} based on MDI continuum observations. Assuming 
the hot wall model (Spruit 1976) and a simple flux tube geom- 
etry, the authors demonstrated this to imply the average size of 
flux tubes increases with {Bi) ///. 

Both continuum and line core intensity specific contrast in- 
creases strongly with decreasing (Z?/) /fi up to a maximum (at 
around 200 and 50 G respectively) before gradually declining 
(Figs. 14 and 15 1. The decline at low {Bi) Ifi aside, the observed 



trends suggests that brightness or temperature excess in both 
the lower and middle photosphere in flux tubes decreases with 
increasing magnetic filling factor (the continuum and line core 
formation height of the Fe I 6173 A line is about 20 and 300 km 
respectively, Norton et al.|2006] l. This may partly be due to the 
increasing average size of flux tubes; larger magnetic elements 
appear darker as lateral heating is less efficient (|Spruit||l976; 



more severe ( |Kobel et al.|2012|l. Th e total radiation flux derived 
from MHD simulations by |V ogler ( 2005 ) exhibits a similar be- 
haviour. This result also agrees with the observations that net- 
work elements appear hotter than facular elements in the lower 



and middle photosphere (Solanki & Stenflo 1984 



Keller et al.|1990| > and G-band bright points are brighter in the 



Solanki 1986 



quiet Sun than in active regions (Roma no et al.|2012) l. 

The decline in both continuum and line core specific con- 
trast at low {Bi) /fi suggests a diminishing heating efficiency of 
the smallest magnetic elements but is more likely due to the in- 
fluence of intergranular lanes on apparent contrast. (As with the 
low contrast towards (Z?/) /jj = G in the contrast versus (Z?/) /fi 
profiles discussed in Sect. |3. 1 1 ) The probable impact of the trun- 
cated coverage of disc positions at low (Z?/) ///on the quality of 
the fits to measured contrast here might have also played a role. 



4. Discussion 

4.1. Comparison with \Ortiz eTai]f2002l 

A comparison of the results reported here for continuum inten- 
sity with those from the similar study based on full-disc MDI 
observations by Ortiz et al. (2002) reveal several notable differ- 



The contrast reported here is generally higher, by as much as 
a factor of about two. The difference between the two studies 
becomes increasingly pronounced with (Z?/) /fi and distance 
from disc centre. 

Approaching disc centre, contrast appears to level off in our 
study (Fig. |7J. In the earlier work, the contrast declines ap- 
proximately linearly towards fi — 1 and there are also marked 
fluctuations about n ~ 0.95 (Fig. 3, |Ortiz et al.|2002| l. 
Near {Bi) /fx = G contrast is negati ve here (Fig. [9] ) but 



positive in the previous study (Fig. 4, Ortiz et al.|2002 



Spruit & Zwaan|1981 Grossman n-Doerfh et al.|1994 1 and mag- 
netic suppression of surrounding convective energy transport is 



- The specific contrast (a proxy of intrinsic contrast given by 
the quotient of contrast at the maximum point on the CLV 
profile and {B/) /fi) presented here ascends with (Z?/) /fi up to 
{Bi) l\i ~ 200 G before descending monotonically thereafter 
(Fig.[l4j>. Ortiz et aL| ( |2002| ) found specific contrast to decline 
approximately linearly with {Bj) //j. (Fig. 8 in their paper). 

- For {Bi) In > 200 G, the specific contrast reported here is 
also nearly double that in the earlier work. 

The lower contrast and specific contrast (for {Bi) jfi > 200 G), 
and difference in CLV towards disc centre reported by Ortiz et al. 
(2002) is likely, as we will show shortly, to be primarily due to 
the misidentification of the magnetic signal adjacent to sunspots 
and pores as network and faculae (discussed in Sect. |2.2| i by 
those authors. Care was taken here to minimise such misidentifi- 
cation. The negative contrast towards {B{) jfi — G found here, 
as argued in Sect. [XT] arises from the resolution of intergranular 
lanes at HMFs spatial resolution. We demonstrate below that this 
difference in spatial resolution also contributes to the opposite 
{Bi) ///-dependence of specific contrast below {Bi) jfi ~ 200 G 
reported here and by Ortiz et al. ( 2002| ). 

To recreate the conditions of the study by |Ortiz et al.| ([2002 ), 
we recomputed continuum intensity contrast and specific con- 
trast from the HMI data set employed here without applying the 
magnetic extension removal procedure, binning the data set spa- 
tially by 4 x 4 pixels to be consistent with MDI's spatial reso- 
lution and (very approximately) transforming measured contrast 
to the corresponding value at MDI's operating wavelength, 6768 
A. This last transformation was carried out by describing quiet 
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Sun, network and faculae as black bodies, and taking an effec- 
tive temperature of 5800 K for the quiet Sun. This allowed us 
to crudely convert contrast measured at 6173 A by HMI into the 
corresponding contrast at 6768 A. It should be noted that this 
is a first-order approximation of the wavelength dependence of 
contrast which ignores the variation of the continuum formation 



dSc 

1999||Norton et al.|2~006) . In Fig.|16|we depictlhe contrast CLV 
profile of network patterns (50 G < (B/) /// < 100 G) and active 
region faculae (500 G < (B/) //j < 640 G) after the application of 
the above procedure. Also plotted are the similarly treated con- 
trast versus (B/) /fj. profiles about disc centre (0.94 < // < 1.00) 
and near limb (0.36 < jj. < 0.50). The specific contrasts from 



height with wavelength dSolanki & Unruh| 19981 ISiitterlin et al. 



this process are illustrated in Fig. 17 
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Fig. 17. Continuum intensity specific contrast as a function of 
{Bf) //i from Fig. 14 (open circles). The cyan, blue and green 



filled circles represent the same quantity obtained by omitting 
the magnetic extension removal procedure, spatially binning the 
data set by 4x4 pixels and converting measured contrast to 6768 
A respectively. The red filled circles represent the result of taking 
into account all these three considerations. The points in each 
series are joined by straight lines to aid the eye. 



HMI's spatial resolution plays a role in the observed decline in 
specific contrast from (B/) /fi ~ 200 G towards {B{) /fj. — G. 

The largest effect is produced by the removal of magnetic 
signal adjoining to sunspots and pores. Hence it is worth con- 
sidering this process more closely. The procedure by which we 
removed these signals inevitably discriminates against active re- 
gion faculae, and various authors have noted lower contrast in 
active region faculae compared to quiet Sun network at similar 
magnetogram signal levels ( |Lawrence et al.p993| Kobel et al. 



201 1 ), but we do not reckon this to be the reason for the greater 
contrast reported here. In spite of the severe steps taken to ex- 
clude magnetic signals adjoined to sunspots and pores, there re- 
mains a fair representation of active region faculae in the mea- 
sured contrasts (Fig. [5 1. Also, at HMI's spatial resolution, net- 
work patterns are largely confined to the lower half of the {Bf) In 
range considered (Fig. 6]) while the difference between the con- 
trast reported here and by |Ortiz et al.| (|2002) is more pronounced 
at higher (Bi) /fi. 



Following the example of earlier studies (e.g. Chapman 
1980), Ortiz et al. ( 2002} fit quadratic polynomials to contin 



uum intensity contrast CLV profiles derived in their study. Here 
we fit cubic polynomials instead because of the different CLV 
near disc centre which we just demonstrated to arise from our 
treatment of magnetic signal adjoined to sunspots and pores. 

We surmise that the fluctuations about disc centre in the con- 
tinuum intensity CLV profiles reported by Ortiz et al. ( 2002| l is 
due to the non-homogeneous distribution of active regions over 
the solar disc. Each unit /i represents a greater radial distance 
towards the disc centre rendering these fluctuations increasingly 
abrupt. Not accounting for magnetic signal adjoined to sunspots 
and pores probably accentuated these fluctuations. 

4.2. Continuum intensity contrast about disc centre 

In this section we compare the magnetogram signal dependence 
of continuum intensity contrast near disc centre reported here 
with that from other recent studies, summarized in Table[2] Apart 
from the work of |Ortiz et al.| ( |200"2") , the CLV of continuum in- 
tensity contrast has been examined in detail previously bypbpka 
eTal] f(fl99"2"l \lWfj and |Lawrence et al.| ( fT$93] >. While we and 
Ortiz et al. ( |2002| l employed full-disc observations made at a 



Excluding the procedure to remove magnetic signal adjoined 
to sunspots and pores produced the largest effect. It produced an 

overall drop in contrast (cyan series, Fig. 16 1 and also repro- from those presented here and by lOrtiz et al 



single wavelength and spatial resolution, these authors collated 
telescope scans made at multiple wavelengths and resolutions. 
The magnetogram signal intervals represented by the contrast 
CLV profiles presented in these papers also differ considerably 

— d2002l). For these 



duced the linear CLV near disc centre observed by |Ortiz et aT] 
p002") (640 G < (B,) In < 800 G panel). The effects of the resiz- 
ing of the data set (blue series) and the projection of measured 
contrasts to MDI's operating wavelength (green series) are rela- 
tively minor. The contrast profiles derived from the application 
of all three procedures (red series) are, in terms of both form 
and magnitude, in general agreement with the profiles covering 
similar disc positions and magnetogram signal levels reported by 
Ortiz et al. ( 2002 1 (Figs. 3 and 4 in their paper). 

Similarly for specific contrast, the result of applying all three 
processes (red series, Fig. [17) resembles the measurements pre- 
sented by Ortiz et al. (2002) (Fig. 8 in their paper) and the great- 
est effect on magnitude came from the omission of the magnetic 
extension removal procedure. Binning the data set down to a 
MDI-like resolution also produced a significant increase in spe- 
cific contrast at low {Bj) jpi. After combining with the other two 
processes, this rendered specific contrast approximately level be- 
low 200 G. The influence of granulation on apparent contrast at 



reasons it is not straightforward to make any quantitative com- 
parisons with the contrast CLV reported by these studies, there- 
fore the focus on results for near disc centre here. 

The negative contrast at low magnetogram signal levels 
found here and by the majority of the compared works coun- 
ters expectations from thin flux tube models, which predict in- 



trinsically bright magnetic features in this regime (Knolker et al. 
1988|). As d iscussed in Sect. |3.1| this is attributed here (and in 



Title et al.|[T992l |Topka et al.||1992] |Kobel et al.|[20lT) to the 



influence of intergranular lanes, an assertion supported by var- 
ious models ( |Title & BergerfT996l |Schnerr & Spruit||201 1) . Of 



the studies compared, only that by Ortiz et al.| ( |2002| l noted pos- 
itive contrasts. Binning the data set here to a MDI-like resolu- 
tion raised the disc centre end of the contrast CLV profile in the 
50 G < {Bi) /fi < 100 G interval into positive territory (Fig. 
[16), suggesting that the positive contrasts reported by |Ortiz et al. 



( |2002) arose from utilising data at a resolution (the lowest of the 
studies compared) where granulation is largely unresolved. 
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Fr om SVST sca ns of active regions, |Title e t al. ( 1992), Topka 

eTaTl ( [1992] [1997) and |Lawrence et al.| ( |1993| ) foun d contrast 
to decline monotonically with magnetogram signal. Lawrence 
et al.| ( |1993[ ) also noted the same with active region data ac- 
quired at the San Fernando Observatory (SFO). |Kobel et al. 
( 201 1) demonstrated the comparatively poorer resolution and 
straylight from pores to be likely culpable for the monotonic de- 
cline observed in the SVST studies. In all the other compared 
works, contrast exhibits a peak, the position of which varied 
from (B,) In ~ 10 G ( |Ortiz et al.||2002) to (B,) ~ 700 G 



(Kobel et al. 201 1), except for the quiet Sun SST scan exam- 
ined by |Schnerr & Spruit ( |2011) , where contrast saturated at 
(Bi) ~ 800 G. Generally, the finer the spatial resolution, the 
higher the position of the peak. Most of these studies were based 
on observations made at the Fe I 6302 A line. Even if we dis- 
count the studies made at other wavelengths, this broad pattern 
is still apparent, ruling out differences in wavelength as the ma- 
jor driver. Ko bel et aL[ ( |201 1) reported similar contrast peak po- 
sitions for quiet Sun and active region scans. Therefore the dif- 
ference between active region and quiet Sun contrast had likely 
little role in the spread in reported peak positions amongst the 
compared works. As shown in the 0.94 < /j < 1.00 panel of Fig. 
16 both resizing the HMI data set to a MDI-like spatial reso- 
lution and omitting the magnetic extension removal procedure 
shifted the position of the contrast versus {Bi) /fi profile maxi- 
mum towards the origin. 

The above comparison points to differences in spatial reso- 
lution and treatment of magnetic signal near sunspots and pores 
as the dominant factors behind the spread in reported depen- 
dence of contrast on magnetogram signal. Indeed, the recent 
MHD simulation work of Rohrbein et al. ( 201 1[ ) suggests that 
the contrast peak at intermediate magnetogram signal levels seen 
in direct observations, but not observed in MHD simulations, is 
a product of the limited spatial resolution. The bulk of the stud- 
ies listed in Table [2] were based on higher spatial resolution data 
than utilised here. For this reason, a more quantitative compari- 
son like we did in Sect. 4. 1 with the findings of Ortiz et al. ( 2002) 
is not workable here. 



The recent works of Berger et al. ( 2007 1 and Viticchie et al. 
(2010), examining the contrast of bright points in the G-band, 
bear tenuous relevance to our study. These studies utilised ob- 
servations made in a molecular band (i.e., neither continuum 
nor line core) where the contrast of magnetic features is en- 
hanced (due to CH depletion, Steiner et al. 2001 Schiissler et al. 
2003 ). More importantly, while this study and the works cited 
in Table [2] considered each pixel a separate entity, Berger et al. 
( |2007) and | Viticchie et al.| ( |2010) examined the overall contrast 
of each individual bright structure. As the body of magnetic fea- 
tures isolated by both approaches differ, the results are not di- 
rectly comparable. Berg er et al.| d2007 ) did however also report 
a pixel-by-pixel consideration of G-band contrast. Scans at four 
disc positions were surveyed. Barring one that appeared anoma- 
lous, the contrast versus magnetogram signal profiles from each 
scan bear general resemblance to ours in terms of form. Notably, 
the profile from the disc centre scan, with a spatial resolution 
of 0.15 arcsec, exhibits a peak at (Bi) ~ 700 G, consistent with 
the broad pattern between spatial resolution and peak position 
described here. 



4.3. Line core intensity contrast 

As pointed out in Sect. |3.1[ the most notable difference between 
the continuum and line core intensity contrasts present here is 



the converse CLV. While continuum intensity contrast is weak- 
est near disc centre and strengthens towards the limb, line core 
intensity contrast is strongest at disc centre and declines to- 
wards the limb. The divergent CLV exhibited by the two sets of 
measurements stem from their rather different physical sources. 
Continuum intensity is enhanced largely in the hot walls of mag- 
netic elements, thus the centre-to-limb increase. The line core is 
formed in the middle photosphere, which is heated either by ra- 
diation from deeper layers (Knol ker et al.|199"T) , or by mechan- 
ical and Ohmic dissipations (Moll et al. 2012). Also stated in 
Sect " 
shape 



3.1 



line core intensity is modulated by line strength and 
and continuum intensity. Excluding variation related to 
continuum excess, line core intensity enhancement arises from 
the influence of the temperature excess in the middle photo- 
sphere and Zeeman splitting on line strength and shape. Given 
the relatively narrow, vertical geometry of flux tubes, as mag- 
netic elements rotate from disc centre to limb, line-of-sight rays 
go from being largely confined to single flux tubes to increas- 
ingly passing into and out of multiple flux tubes, especially in 
densely packed facular regions ( jBiinte et al.|1993 I. Line core in- 
tensity contrast decreases towards the limb as the contribution 
to the spectral line from magnetic elements diminishes from ab- 
sorption in the non-magnetic part of the atmosphere transversed 
by the rays (Solan kTet al.|1998 i. Another probable cause of the 
centre-to-limb decline is the spatial displacement of the line core 
with respect to the corresponding continuum towards the limb 
caused by the difference in formation height and oblique viewing 
geometry (Stellmac her & Wiehr|199l" 2001 1. The line core in- 
tensity enhancement arising from temperature excess in the mid- 
dle photosphere and Zeeman splitting discussed here is not to be 
confused with that from the centre -to-lim b broadening of the Fe 

which arises from the 



2.2.2 



I 6173 A line mentioned in Sect, 
viewing geometry independent of magnetic field. 

We recomputed line core intensity contrast, this time nor- 
malizing the line core intensity images by the corresponding 
continuum intensity images prior to data reduction. The result, 
the line core residual intensity contrast, is essentially the com- 
ponent of line core intensity contrast arising from line weaken- 
ing in magnetic features alone. Line core intensity and resid- 
ual intensity contrast values can be compared directly. (Line 
core residual intensity contrast equates to the line core inten- 
sity contrast we would get from scaling the line core intensity 
of just the network and faculae pixels by (/c,qs) lie ) The line 
core residual intensity contrast CLV profile of network patterns 
(50 G < (Bi) In < 100 G) and active region faculae (500 G < 
(Bi) IfJ. < 640 G), and contrast versus {Bj) /fi profile about disc 
centre (0.94 < fi < 1.00) and near limb (0.36 < fi < 0.50) so 
derived are plotted along with the corresponding line core inten- 
sity contrast profiles from Figs. [8] and 10 in Fig. 18 Line core 



residual intensity specific contrast is plotted along with the line 



core intensity specific contrast from Fig. 15 in Fig. 19 



As expected, excluding the contribution of continuum inten- 
sity enhancement to line core intensity contrast results in a more 
pronounced centre-to-limb decline (50 G < < 100 G 

and 500 G < (B,) //j < 640 G panels, Fig. [T8). This decline 
in the line core residual intensity contrast CLV profiles is also 
consistent with the spectroscopic observations of Stellmacher & 
|Wiehr| (1979) and |Hirzberger & Wiehr| ( |2005) . Comparing line 



core intensity and residual intensity contrast (Figs. 18 and 19 1 
the broad similarity, especially near disc centre, imply that line 
core intensity contrast is dominated by the contribution from line 
weakening. Even towards the limb, where magnetic line weaken- 
ing is at its weakest, it still comprises a significant proportion of 
observed line core intensity contrast. The potential implications 
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Fig. 19. Line core intensity (open circles) and corresponding 
residual intensity (filled circles) specific contrast as a function 
of (Bi) //u. The circles are joined by straight lines to aid the eye. 



of this result and the opposite CLV exhibited by continuum and 
line core intensity contrast on facular contribution to solar irra- 
diance variations will be discussed in Sect. 14.41 

Very few studies examining the disc position and magne- 
togram signal dependence of line core intensity contrast exist 
in the literature. |Frazier| ( |19711 l, [Lawrence et al.| (fl991) and |Title| 
et al. ( 1992 1 measured intensity contrast at disc centre in the core 
of the Fe I 5250 A, Fe I 6302 A and Ni I 6768 A lines, respec- 
tively. Results from different lines are not directly comparable 
but it is still encouraging to see that in terms of magnitude and 
general trend with magnetogram signal, the results from these 
earlier studies express broad agreement with ours. The only no- 
table exception is the steep monotonic decline at {Bf) > 600 G 



reported by |Title et a l. ( 1992), which largely persisted even after 
the authors masked out pores. The decline coincides with a simi- 
lar trend in intensity contrast measured at the nearby continuum, 
suggesting a greater relative influence of the continuum excess 
to magnetic line weak ening at the Ni I 6768 A line compared 
to the Fe I 6173 A line. |Frazier| ( |l97l] l also made measurements 
away from disc centre, but due to data scatter the author could do 
no better than to present a schematic representation of the CLV 
through visually fitted linear functions . | Walton| ( | 1 987| l reported 
measured line core residual intensities for eight magnetically 
insensitive lines at various disc positions. There were however 
only a small number of scattered measurements for each line, 
and a relatively narrow range of disc positions (/i > 0.6). To our 
knowledge, the present study is the first to examine line core in- 
tensity contrast employing full-disc observations. Together with 
the relatively fine resolution and low noise level, this allowed us 
to describe the magnetogram signal dependence and especially 
the CLV with greater accuracy and detail than in the previous 
efforts. 



4.4. Facular contribution to variation in solar irradiance 

Our finding that intensity contrast in the line core is dominated 
by line weakening rather than continuum intensity enhancement 
and exhibits the opposite CLV as continuum intensity contrast 
has potential implications on facular contribution to total solar 
irradiance variation. Total solar irradiance (TSI) variation is the 
sum manifest of variation in the continuum and spectral lines. 



The converse CLV of the continuum and line core intensity con- 
trasts reported here imply a different time variation of the con- 
tribution by the continuum and spectral lines to solar irradiance 
as magnetic features rotate across the solar disc. 

In Fig. 20 1 we show the level 2.0 hourly TSI measure- 
ments from the DIARAD radiometer (Dewi tte"et al.|20 04 ) on the 
SoHO/ VIRGO instrument ( [Frohlich et al.|1995| l for the 22-day 
period of August 19 to September 9, 1996 (open circles). Active 
region NOAA 7986 rotated into view on August 23 and out of 
view on September 5, and was the only active region on the solar 
disc for the duration of its passage across the solar disc. The pe- 
riod is otherwise relatively quiet. This data therefore allows us to 
chart variations in TSI arising mainly from the passage of a sin- 
gle active region across the solar disc. The dip around August 29 
corresponds to active region NOAA 7986 crossing disc centre. 
Even with darkening from sunspots present in the active region, 
the nadir of the dip is ~ 0.2 WirT 2 above the level before August 
23 and after September 5. This suggests an overall positive con- 
tribution to TSI variation by the faculae in NOAA 7986 when it 
was near disc centre (first pointed out by Fligge et al. 2000). 

A schematic representation of facular contribution to varia- 
tion in TSI during the passage of NOAA 7986 across the solar 
disc, depicted in Fig. 20:, was derived as follows. 



The DIARAD TSI data was interpolated at 0. 1 day intervals 
and the result smoothed (via binomial smoothing, Ma rchand| 



& Marmet 1983 1. Sunspot darkening was estimated from 



the Photometric Sunspot Index (PSI) by Chapman et al. 
( 1994) based on full-disc photometric images acquired with 
the Cartesian Full Disk Telescope 1 (CFDT1) at SFO A 
quadratic polynomial was fitted to the PSI values from the 
period of interest excluding the points where PSI = (i.e., 
no sunspots in view). Taking a value of 1365.4 WirT 2 for the 
total irradiance of the quiet Sun, the fit was converted from 
units of parts per million to WirT 2 and subtracted from the 
DIARAD data. This value for the total irradiance of the quiet 
Sun is given by the average TSI at the last three solar min- 
ima stated in version d41_6 2_1204 (dated A pril 2, 2012) of 
the PMOD TSI composite ( |Frohlich|[2000| . The DIARAD 
data, after this treatment, represents variation in TSI largely 
from faculae in NOAA 7986 alone. In Fig. 20 1 we plot the 



DIARAD data after interpolation and smoothing (red dot- 
ted curve), and after removing sunspot darkening (red solid 
curve) along the original measurements. In Fig. 20 3 we show 
the PSI (open circles) and the quadratic polynomial fit to the 
non-zero segment (curve). 

The trajectory of NOAA 7986 during its passage across the 
solar disc was estimated from 142 level 1.8 5-minute MDI 
magnetograms ( Sch errer et al.|1995] ) on which the active re- 
gion was entirely in view (i.e. not only partially rotated into, 
or partially rotated off the solar disc). Taking the unsigned 
magnetogram signal, the magnetograms were binned spa- 
tially by 16 x 16 pixels. For each binned magnetogram, the 
position of NOAA 7986, in terms of fi, was estimated from 
the mean position of the five pixels within the active region 
with the strongest signal. The trajectory of NOAA 7986 is 
then given by the quadratic polynomial fit to these estimates. 
In Fig. [20fc we show the estimated position of NOAA 7986 
in the magnetograms (open circles) and the quadratic poly- 
nomial fit (curve). 

Facular contribution to variation in TSI was very approxi- 
mately modelled from the empirical relationships describing 
contrast as a function of fi and (B/) //j derived in this study. 
Assuming a power law distribution of (B/) //j. with a scal- 
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ing exponent of -1.85 (Parnell et al. |2009[ ), we evaluated 

Eqs. 3 and 4 scaled by (^/15)~ ' and integrated over 
{Bi) fjx = 15G to 800 G, at 0.1 day intervals taking (i from 
the trajectory of NOAA 7986 estimated earlier. The result- 
ing time series were then scaled by /i (to correct for the CLV 
of projected area on the solar disc) and the limb darkening 
function from Fouk al et al.| (|2004). Given the approximate 
manner of this derivation, the results are non-quantitative. 
However for this analysis it is not the actual values but the 
temporal trends that is important. Here we had approximated 
the active region as a point object and so cannot include vari- 
ation arising from the active region being only partially vis- 
ible as it rotates into and off the solar disc. This derivation 
is therefore only valid, and confined to, the period where 
NOAA 7986 was entirely in view in MDI magnetograms. 



The treated DIARAD data, giving TSI variation largely from 
faculae in NOAA 7986 alone, is plotted in Fig. 201 (red curve) 
along the conjectures based on the observed intensity contrast 
in the continuum (blue dashed curve) and line core (blue dotted 
curve). To compare how they varied with time, we subtracted the 
mean from and normalized each time series by the area bounded 
by the curve and the zero level. In Fig. 



20 5 we show the multiple 



linear regression fit of the continuum and line core models to the 
DIARAD series (blue solid curve), which showed a much better 
agreement to it than either model. Facular contribution to varia- 
tion in solar irradiance appears to be strongly driven by intensity 
contrast in both the continuum, and spectral lines, which derives 
largely from line weakening in magnetic elements. This com- 
plies with the observation that solar irradiance variations over 
the solar cycle seems to be significantly influenced by changes 
in spectral lines ( Mitcheff & Livingston| 1 99 1 [ |Unruh et"aT 1 999 



Preminger et al.|2002) > 



Various studies examining the photometric contrast of facu- 
lae, as identified in calcium line images, reported positive val- 
ues near disc centre ( |Lawrence et al.|1 988 ; Wal ton et aH2 003 1 
Foukal et al. 2004). The apparent divergence between the find- 
ings of these works, and the largely negative contrast near disc 
centre reported here (Fig. [9]) and in the similar studies listed in 
Table [2] where network and faculae were characterized by the 
magnetogram signal, arises from the different selection methods 
( |Ermolli et al.||2007] ). Significantly, the results of |Foukal et al. 
(2004 1 were based on near total light broadband measurements. 
From the statistical analysis of an extensive catalogue of full- 
disc photometric observations obtained with the SFO/CFDT1, 



Walton et al.| (|2003[) found contrast, in a continuum passband 



(i.e., negligible spectral line contribution) in the visible red, av- 
eraged over each facular region (as identified in the Ca II K) 
to be positive (~ 0.005) near disc centre. The authors took this 
result to indicate the distribution of flux tube sizes in facular re- 
gions is biased to the low end (i.e., more small bright flux tubes 
than large dark ones), which is reasonable considering the typ- 
ical distribution of magnetogram signal (i.e., much more weak 
signals than strong, Parnell et al.|20 09 ) and the observation that 
the average flux tube size is greater where magnetogram sig- 
nal is greater (Ortiz et al. 2002). This suggests the continuum 
component of facular contribution to TSI variation is on aver- 
age positive at disc centre in the visible red. From this and the 
analysis presented above, we conclude that the apparent overall 
positive contribution to TSI variation by the faculae in NOAA 
7986 when it was near disc centre is, at least in part, due to the 
prevalence of smaller flux tubes and line weakening. 



5. Conclusion 

Here we have presented measured network and facular intensity 
contrast in the continuum and in the core of the Fe I 6173 A 
line from SDO/HMI full-disc observations. We studied the de- 
pendence of the contrast on disc position and magnetogram sig- 
nal, represented by fi and {Bi) /fj. respectively. Specifically, we 
derived empirical relationships describing contrast as a function 
of fi and {B/} /ju, and specific contrast (contrast per unit {B/} //i, 
representative of intrinsic contrast) as a function of (B/) //j. 
This study exploits the unprecedented opportunity offered by 
the SDO mission to examine co-temporal full-disc observations 
of magnetic field and intensity at a constant intermediate spa- 
tial resolution (1 arcsec), relatively low noise, and without at- 
mospheric interference. The quality of the data allowed us to 
examine intensity contrast for a larger sample at greater accu- 
racy and detail than previous, similar studies, especially in the 
case of the line core. These results constitute stringent obser- 
vational constraints on the variation of network and facular in- 
tensity with disc position and magnetogram signal in the low 
and middle photosphere. By constraining atmospheric models 
of network and faculae, these results should be of utility to solar 
irradiance reconstructions, especially as HMI data will increas- 
ingly be used for this purpose. Given this is the first study of its 
kind to examine the entire solar disc in both the continuum and 
line core, it should also be useful to reproduce these results in 
models of magnetic flux concentrations. 

There are significant discrepancies in the continuum inten- 
sity contrast reported here and from earlier studies. In this study 
we had taken steps to account for magnetic signal in the periph- 
ery of sunspots and pores, arising from their magnetic canopies 
and the influence of straylight, which can easily be misidentified 
as network and faculae. From a comparison with the findings 
of past efforts, including a recomputation of the results obtained 



here recreating the conditions of the similar study by Ortiz et al. 



(2002), we showed differences in resolution, and treatment of 
magnetic signal adjacent to sunspots and pores to be the likely 
main factors behind the spread in reported results. The apparent 
radiant behaviour of network and faculae elements is strongly 
modulated by spatial resolution (|Title & Bergerfl 996; Rohrbein 
|et al.|20lTj |Schnerr & Spruit|2011[ ). An understanding of its in- 
fluence is necessary for the proper interpretation of direct mea- 
surements. 

In terms of magnitude, trend with magnetogram signal and 
in particular the CLV, the results obtained here in the contin- 
uum and line core differ considerably. While continuum inten- 
sity contrast broadly ascends towards the limb, line core inten- 
sity contrast is greatest near disc centre and diminishes from disc 
centre to limb. The divergence between both sets of measure- 
ments arises dominantly from spectral line changes due to heat- 
ing in the middle photosphere and Zeeman splitting in magnetic 
features, and the different mechanisms by which apparent con- 
trast vary with viewing geometry going from disc centre to limb. 
From a simple model based on the empirical relationships be- 
tween contrast, and yiz and (Bi) jyt derived here we confirmed that 
facular contribution to variation in solar irradiance is strongly 
driven by both continuum excess and spectral line changes. 

The specific contrast in both the continuum and line core ex- 
hibit a marked decline with increasing magnetogram signal, con- 
firming that network elements are, per unit magnetic flux, hotter 
and brighter than active region faculae. The different radiant be- 
haviour of network and faculae, not accounted in present mod- 
els of solar irradiance variations, would be an important factor 
to consider for more realistic modelling. This observation also 
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implies that secular changes in solar irradiance may be consider- 
ably larger than what some models of solar irradiance variations 
have suggested, given the variation in the number of small-scale 
magnetic elements on the solar disc is a prime candidate driver 
of secular changes (Solanki et al. 2002). For example, the model 



employed by Krivova et al. ( 2007 i to reconstruct variation in TSI 
from 1700 and by Vieira et al.|( |201 l|l for over th e Holocene as- 
sumed the faculae contrast model by |Unruh etaL ( 1999| > for both 
network and faculae. This renders network and faculae with sim- 
ilar magnetic filling factors equally bright and therefore possibly 
underestimate the contribution by network to secular variation. 
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Fig. 6. 100 x 100 arcsec insets of a quiet region near disc centre (p. > 0.99, left) and active region NOAA 11187 (0.82 < fi < 0.91, 
right) of the continuum intensity image from April 15, 201 1. The cyan contours indicate the boundary of network and faculae. The 
blue, green and red contours correspond to (Bi) /fi = 100 G, 280 G and 500 G, respectively. 
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Fig. 7. CLV of continuum intensity contrast, Cj c over eight (B/) j\i intervals. The filled circles and error bars represent the mean and 
standard deviation of C/ c grouped by yU in bins 0.05 wide. The red curves are third-order polynomial fits to the filled circles and the 
blue curves are the cross-sections of the surface fit to C/ c at the mean {Bf) jyi within each interval (largely hidden due to the close 
agreement), extrapolated to fi = (dotted segments). The horizontal dashed lines denote the mean quiet Sun level. 
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Table 2. Continuum intensity contrast dependence on magnetogram signal near disc centre obtained in the present and from earlier 
studies. 



Reference 


Instrument 


Target* 


Resolution [Arcsec] 


Wavelength [A] 


Results 


Present 


SDO/HMI 


H > 0.94 c 


1 


6173 


Negative near (B,) /fi = QG, Peak at (B,) In ~ 200 G. 


1,2,3,4 


SVST 


AR 


>0.3 


Various'' 


Negative near (B t ) = G. Monotonic decline with (Bi). 


3 


SVST 


QS 


>0.3 


6302 


Negative near (B,) = G. Peak at (B,) ~ 400 G. 




SFO/SHG" 


QS 


> 1 


6302 


Negative near <fi ; > = G. Peak at 100 G < (B,) < 200 G. 






AR 






Negative near (Bi) = G. Monotonic decline with {Bi). 


5 


SoHO/MDI 


li > 0.92 c 


4 


6768 


Positive near (Bi) /ji = G. Peak at (B,) In ~ 100 G. 


6 


SST 


AR 


0.15 


6302 


Negative near (B,) = G. Peak at (B t ) ~ 650 G. 


7 


Hinode/SOT 


QS 


0.3 


6302 


Negative near (B,) = G. Peak at (B,) ~ 700 G. 






AR 






Negative near (B,) = G. Peak at (B,) ~ 700 G. 


8 


SST 


QS 


0.15 


6302 


Negative near (B,) = G. Saturate at (B,) ~ 800 G. 




Hinode/SOT 


QS 


0.3 


6302 


Negative near (B,) = G. Peak at (B,) ~ 500 G. 



References. (1) Title et al. ( 1992); (2) Topka et al. 


( 1992); (3) Lawrence et al. ( 1993) 


|& Scharmer|(|2010f; (7)|Kobel et al.|l|201 \); (8)|Sc 


hnerr& Spruit|p01 If. 



Notes. (a) San Fernando Observatory (SFO) 28-cm vacuum telescope and vacuum spectroheliograph (SHG). ib) AR and QS denote active region 
and quiet Sun respectively. (c) Segment of full-disc observations considered. (rf) 5250A, 5576A, 6302A and 6768A. 
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Fig. 9. Continuum intensity contrast, C/ c as a function of (Bi) //j. over eight /u intervals. The filled circles and error bars represent the 
mean and standard deviation of C/ c grouped by {Bi) /ft in 40 G bins. The red and blue curves and the horizontal dashed lines have 
the same meanings as in Fig. [7] 
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Fig. 10. Similar to Fig.|9j but for line core intensity contrast. 
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Fig. 16. Selected continuum intensity contrast CLV (left) and contrast versus (B/) /fi (right) profiles from Figs. |7]and|9](open circles 
and dotted curves). The selected profiles correspond to quiet Sun network (top left), active region faculae (bottom left), disc centre 
(top right) and near limb (bottom right). The cyan, blue and green series denote the profiles obtained by omitting the magnetic 
extension removal procedure, spatially binning the data set by 4 x 4 pixels and converting measured contrast to 6768 A respectively. 
The red series indicates the results of taking into account all these three considerations. The circles represent the mean of measured 
contrast binned as in the referenced figures and the curves the corresponding third-order polynomial fits. The dashed lines mark the 
mean quiet Sun level. 
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Fig. 18. Selected line core intensity (open circles and dotted curves) and corresponding residual intensity (filled circles and solid 
curves) contrast CLV (left) and contrast versus (Bi) /fi (right) profiles. The circles represent the mean of measured contrast binned 
as in Figs.[8]and 10 and the curves the corresponding cubic polynomial fits. The selected profiles correspond to quiet Sun network 
(top left), active region faculae (bottom left), disc centre (top right) and near limb (bottom right). 
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Fig. 20. a) Total solar irradiance, TSI from DIARAD on SoHO/VIRGO (open circles) for the period of August 19 to September 
9, 1996. The dotted and solid red curves represent the interpolated and smoothed version before and after subtracting sunspot 
darkening, b) Photometric Sunspot Index, PSI (open circles) by Chapman et al. ( 1994) and the quadratic polynomial fit to the non- 
zero points (curve). The dashed line marks the zero level, c) Position of NOAA 7986 as estimated from MDI magnetograms (open 
circles) and the quadratic polynomial fit (curve), d) TSI minus sunspot darkening (red curve) and model of facular contribution 
to TSI based on observed intensity contrast in the continuum (blue dashed curve) and line core (blue dotted curve), each mean 
subtracted and normalized by the area bound between the time series and the zero level (black dashed line), e) TSI minus sunspot 
darkening (red curve) and the multiple linear regression fit of the continuum and line core models (blue solid curve). The dotted 
lines running down all the panels mark the period where NOAA 7986 was entirely on the solar disc in MDI magnetograms. 
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